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ABSTRACT 

This  paper  is  concerned  with  the  application  of  the  Ritz-Galerkin 

method  to  the  numerical  solution  of  singular  boundary  value  problems  of  the 

type  arising  when  Poisson's  equation  on  a domain  with  cylindrical  or  spherical 

symmetry  is  reduced  to  a one-dimensional  problem.  The  objective  is  to  derive 

a priori  L^  and  L -norm  estimates  for  the  error.  The  difficulty  is  that 
2 

these  norms  are  not  natural  norms  for  the  reduced  problem.  With  the  aid  of 
B-splines  we  prove  some  nonstandard  approximation  - theoretic  results  and 
use  these  to  derive  the  desired  error  estimates.  Some  numerical  results 
are  presented . 

AMS  (MOS)  Subject  Classifications:  65L10,  65N15,  65N30,  41A25 

Key  Words:  Singular  problems.  Weighted  spline  projections,  Rayleigh- 

Ritz-Galerkin  methods 

Work  Unit  Number  7 (Numerical  Analysis) 

EXPLANATION 

The  Ritz-Galerkin  (finite  element)  method  is  a widely-studied  metho 
for  obtaining  numerical  solutions  to  differential  equations.  Much  is  known 
about  the  convergence  behavior  of  the  approximate  solution  when  all  the 
coefficients  in  the  differential  equation  are  regular.  This  report  investi- 
gates the  application  of  the  method  to  boundary  value  problems  in  one  space 
dimension  with  a singular  coefficient  of  the  type  which  arises  when  symmetry 
permits  the  reduction  of  a partial  differential  equation  to  an  ordinary 
differential  equation.  The  basic  result  is  that  the  singular  coefficient  does 
not  degrade  the  performance  of  the  method.  Some  numerical  examples  are  given. 

f 
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RITZ-GALERKIN  METHODS  FOR  SINGULAR  BOUNDARY  VALUE  PROBLEMS 


Dennis  Jospersen 

51 . Introduction 

When  Poisson's  equation  -fiu  = f is  encountered  on  a domain  with  cylindrical  or 
spherical  symmetry,  one  can,  if  the  data  depend  only  on  the  radial  coordinate,  make  a 
change  of  variables  to  reduce  the  problem  to  a one-dimensional  problem,  albeit  at  the 
expense  of  introducing  a singularity  in  the  equation.  This  report  investigates  the  Ritz- 
Galerkin  process  for  numerically  solving  the  one-dimensional  reduced  problem,  witli  the 
objective  of  deriving  optimal  error  bounds  of  the  type  normally  encountered  in  finite 
element  analysis. 

In  pursuing  Diis  work,  the  motivating  ideas  arc  roughly  the  following.  Assuming 
smooth  data,  the  solution  u to  Poisson's  equation  is  smooth,  hence  is  well-approximable 
by  piecevnse  polynomials.  One  would  thus  expect  a Rit.-’.-Gal  erk  i n approximation  to  u 
from  a i>iecewise  polynomial  space  to  be  "optimally"  close  to  u.  Unfortunately,  the 
singularity  in  the  one-dimensional  problem  frustrates  the  usual  analysis  unless  one 
either  adds  to  the  finite-dimensional  space  functions  which  match  the  b«'havior  of  the 
Green's  function  or  makes  estimates  in  a natural  weighted  Sobolev  norm  (12),  (13]). 

There  have  been  some  recent  works  that  liave  made  progress  on  this  question.  Duixrnt  and 

2 

Wahlbin  (10)  analyzed  the  problcmi  -{a  u')'  + qu  - f on  (0,1)  where  q(x)  ^ q^  > 0 

for  0 £ X £ 1 and  a(x)  c C^[0,1].  Tliey  were  able  to  show  the  Ri  tz-Galcrkin  process 

produces  an  approximate  solution  u.  wliich  is  optimally  close  to  u in  tlie  L_(0,1) 

h 2 

norm.  Their  proof  seems  to  depend  strongly  on  the  conditions  a f C^(0,11  and  q(x) 
strictly  positive.  In  anotlier  recent  work,  de  Iloog  and  Weiss  [7]  have  analyzed  tlie 
application  of  collocation  methods  to  singular  boundary  value  problems  such  as  the  ones 
considered  here.  In  a difficult  piece  of  analysis  they  show  that  collocation  gives, 
roughly  speaking,  the  same  results  for  singular  boundary  value  problems  as  for  nonsingular 


Spon.sorod  by  tlie  United  .Stato.s  Army  under  Contract  No.  DAAG29-75-C-0024  and  the  National 
Science  Foundation  under  Grant  No.  MCr.7S-173R5. 


problems  (for  wliich  see  (6)  or  114)).  Their  results  gave  further  impetus  to  the  hope 
that  the  Ritz-Galerkin  method  would  give  the  same  results  for  singular  problems  as  for 
nonsinqular  problems. 

This  paper  divides  into  three  parts.  In  section  2 we  prove  some  results  of  a 
purely  approximation-theoretic  nature  concerning  the  convergence  of  some  weighted 
projection  operators.  In  section  3 those  results  are  applied  to  give  a priori  error 
estimates  in  and  for  the  basic  problem.  Some  of  ttiese  results  are  extended 

to  nonlinear  problems  in  section  4,  and  some  numerical  examples  are  presented. 

I would  like  to  tliank  Carl  de  Boor  for  several  very  helpful  conversations  in  the 
course  of  this  work. 
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82.  Convergence  of  Weighted  Projections 


Let  0 = X < X < • • • < X “ 1 be  a partition  of  the  interval  I = (0,1).  Define 
U 1 N 


I,  (x.  ,,x.),  h.  :=x.  -X.  ,,  h :=  max  h,.  For  J an  interval,  let  II,  (J)  denote 
1 1-1  11  1 1-1  1 )c 

1 


the  set  of  polynomials  of  order  k (degree  < k)  on  J.  Let  0£v<k-l  be  integers. 


1/et  = sj'  :=  {v  e h'^(I)  : vl  t II,  (I.),  1 < i < n),  where  h'’(I)  denotes  the 
k,v  k 1 — — 


usual  Sobolev  space  of  functions  with  v weak  derivatives  on  1,  equivalently,  the 


V”1  (v“l)  (v) 

space  of  functions  in  C (I)  with  v absolutely  continuous  and  v e L^d). 


Define  the  sequence  t = 


n+k 


= (x-, . . . ,x  ,x  , . . . ,x  , . . . ,x  , . . . ,x  ) where 
U U I 1 N N 


X are  repeated  k times  and  x.  is  repeated  k - v times  for  1 < i < N - 1. 

0 u 1 - - 


It  is  well-known  that  S has  dimension  n = N(k  - v)  + v,  and  a convenient  basis 

,n 


consists  of  the  normalized  B-splines,  {C.},.  We  have  supp  B.  C (t.,t.  .],  B.(x)  > 0, 

1 1 1—1  i+)t  1 — 


.n 


and  B^  (x)  - 1 for  all  x c I.  For  these  and  other  facts  about  B-splines,  see. 


e.g.  |5). 

Let  6 > 0.  Define  a map  Do  - 1*0  . • L (I)  ->  via  / x*^(u  - P.u)';  dx  = 0 


for  aJ  1 f S^,  For  B = 0,  P u is  the  L -projection  of  u onto  and  estimates 

p 2 


of  the  error  u - P^u 

p L 


are  well-known  for  p = 2 and  p = Our  objective  is 


to  derive  a priori  estimates  tor  ||u  - 

p L 


for  6 > 0 and  p = 2," 


Wo  will  make  the  following  assumption  on  the  mesh:  there  are  constants  M > 0, 

y > 1 independent  of  h sucli  that 


(2.1) 


x./x,  < M(i/i)^  for  1 < i < j . 


Define  the  global  mesh  ratio  by  M,  :=  max  (t.  , 

t . . i+k 

= 1.3 


t.)/(t.  , - t.).  We  can  now 
1 D-lk  3 


state  the  main  result  of  this  section. 


Theorem  2.1.  (1)  P„  is  bounded  as  a map  on  L_ , 

p 2 


iPgIl  2 5.  C(k,P,M,Y)  . 


(2)  P„  is  bounded  as  a map  on  L , 
p n> 


< • C(k,B,M,Y) 
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AS  an  immediate  corollary,  wo  have 

Corollary  2.2.  (1)  |!p„u  - u ||  . < C(k,P,M,Y)  inf  ||u  - u.  II  . . 

p n L_ 

2 u.  2 

h 

(2)  ||PgU  - u||  ^ • C(k,6,M,Y)  • inf  ||u  - UjJI  j if  u e L^(I)  . 


The  proof  of  the  Corollary  just  uses  the  fact  that  P is  a projector,  in  addition 

p 

to  Theorem  2.1. 

The  proof  of  Theorem  2.1  will  follow  as  a sequence  of  lemmas.  The  basic  approach 

is  to  write  down  the  matrix  representinq  the  map  P with  respect  to  the  B-spline 

p 

basis  and  then  make  some  careful  estimates  on  Llio  matrix  pioblem. 


% 

Lemma  2.3.  The 

assumjJtion  (2.1) 

impli os 

•1 

(2.2) 

^4k  ^ 

C(M,Y,k)t^  for  k 4 

1 £ i < n , 

f* 

< 

(2.3) 

''j4k'^*^i4k 

■ C(M,Y,k) (j/i)^  for 

1 < i < j s n . 

* 

Proof.  We  have 

^i  *'i(k-v)+j 

for  1 < i < N - 1, 

1 4 V < j £ k,  Xq  = tj 

for 

♦ 

M 

1 < j < k,  and 

^N  " *^N(k-v)4  j 

for  1 4-  V < j < k 4 

V.  lloiice  t.  = X 

1 p(i) 

where 

p ( j ) = N if  n 

4 1 < i £ n 4 k , 

l>(i)=0  if  l<i 

< k,  and  p(i)  = [_( i - 

V - l)/(k 

if  k 4 1 < i < 

n,  where  |_  J 

denotes  tlic  great  i st 

integer  function.  Thus 

t,  /t,  = X ,/x  < M(p(i 

Hk  1 P(r+k)  p(])  - 

4 k)/p(i)  )^. 

For  k 4 1 

i < ri  - k we 

liave 

p(i  t k)  _ - lj/(k  - v)J  ^ . , 

P(i)  |^(i  - V - l)/(k  - v)j  -• 


and  for  n < i t k < n • k we  tiavc 


P(i  t k)  ^ 


‘■p(i)  ~ "[(i  - V - l)7(k  - v)J-  max(l,¥  - k)  ^ ^ ^ ^ • 

Thus  1,  = C(M,Y,k),  which  Jiroves  (2.2). 

Now  consider  (2.3).  If  i = 1,  t.  ,/t.  , < x,/x  < Mj^.  If  i > 1, 

jik  itk  — 3 1 

•"jHk^^i  + k " ’‘p(j+k)^’‘p(isk)  - t k)/p(i  ^ k))^  < M(  (j  - 1 4 k ' v)/(i  - 1))^  < M((l  + k)  j/i)^  , 

which  proves  (2.3)  with  C(M,Y,k)  = M(1  + k)^. 
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Lemma  2.4.  For  1 < i < n , 


o R 

/ X B^(x)dx  > C(H,Y,k,B)t^^j.  / B.  (x)dx  . 


Proof.  For  k + 1 ^ i ^ n,  we  have 


/ X B.  > t.  / B.  = / B.  > C(M,Y,k)-  / B. 


by  (2.2). 


1 f.  i f.  k,  / = / ^ I ^ ^ x^B^.  On  10, x^^]  , B^  is  a polynomial 


of  order  k.  By  the  equivalence  of  all  norms  on  a finite-dimensional  space, 

^1  ^1 
/ X B.  ^ C(k,  B)x^  / B^  . 


’'i  t.  , 

/ x'^P.  > C(k,B)x,  / B.  + xV  / B. 

0 

> min(l,C(k,B))  • xf  / B. 

= C (k,P)  (t.  ^,/t.^  / B. 

k+1  i+k  i+k  •'  1 

> C'(k,B)  C(M,Y,k)"®i"®^t®^.  / B. 

— i+k  •'  1 

> C'(k,B)  C(M,Y,k)'®k'^^t®  , / B. 

~ i+k  ' 1 

= C(M,Y,k,B)t^^j^  / B.  . 

Lemma  2.5.  For  1 < i < n,  there  exists  f.  . e L (1)  with  supp  f.  „ C [t.,t.  .), 
~~  i,B"  i,B— 1 i+k 

lUi^pljj  1 C(M,Y,k,P)  - t^)  and  / x^f^  g(x)B^(x)dx  = for  1 £ j f.  n. 

Proof.  By  M),  v?e  must  consLiiict  f.  „ so  that  x^f.  „(x)  = f^*'^(x)  where 

i , B 1 , B 

(k1 

f 0 k-fold  at  t^ 


(«,(t)  :=  (t  - t.^^)  ...  (t  - V,.,)/(k  - 1)1) 


(2.5)  f = ( 0 at  all  t . f (t. ,t.  , ) 
I J 1 x+k 

I i/i.  k-fold  at  t.  , 

t.  1 i+k 
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If  > 0 we  can  simply  use  the  construction  in  14)  and  get  the  bound 


< t7®D.  (t.^.  - t, ) ^ 

— 1 k i+k  i 


= <w^i>\vV^tk  - 


< C(M,Y,k,3)D^t-®^(t.^^  - t.)-'  . 


if  '"i  ~ ^ slightly  different  construction  suffices.  Let  * :=  t.^j^/4,  and 

let  G be  defined  by 


G (X)  = 


for  X < 5 


for  X > 3iS 


(x  - 6)/26  for  fi  < X 35 


Let  S3  f C (K)  with  supp  'f  £ ss  ^ 0,  / V>  = 1.  For  any  c > 0,  let 

(x)  :=  s5  (x/E)  . 

Now  let  c = 6/2  and  lot  G = G**  * ^ , i.e.,  G(x)  = / G^(y)s5  (x  - y)dy.  Then 

G r c (K)  with  G(x)  = 0 for  x < e,  G(x)  = 1 for  x > — t.  , . 

— — 8 l+k 

Wo  have  |g*^’ (x)  I _<  E ^ ||g*||  ^ 1 Ik  ^ ^ ' II  L 

"I  ■^l  ^ 

f (x)  G (x)  If/ , (x)  . Then  f satisfies  (2.5).  Let  f.  (x)  :=  x~®f^’''(x).  To  estimate 

1 1 # P 

II  f.  „||  , we  note  that  f.  (x)  =0  for  x < E,  hence 

i r P li  1 # P 


le,ial;,c-'  | u,  | |o‘J>  w | 

j=0  J 

< b\-\  y f!)  t^-^8^:>. 


< c(k,6)t:J^(t.^,  - ' 


and  the  proof  is  complete. 


For  u(x)  given,  P^u  t is  obtained  by  solving  the  linear  system 


=oS  ■ S ■ 


kJ 
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where  (G„)  . . : = 

r 8 

j X B.B .dx,  u.  := 

/ x^uB.dx. 

Define 

0 13 

J X 3 1 

D :=  diag((/ 

)'/"),  (Gp)..  ==  j 

-1  -1  -2 

( x^B.B.dx. 
X 3 

Tlicn  1 

equivalent  to  a 

= (D  D)D  u. 

But  we  have 

,1/2 


DGgD  - ^qi  so  G^a  = u is 


iPpuiiL  =115:  Vi^L  ^llsIL  -"d  iII^IIl 


i -. 

?■' 


Hence 

(2.7) 


Iv'Il  ill'^'V'^L  INI 


Lemma  2.6.  There  exists  a constant  = K^(M,Y,k,B)  such  that  for  all  a = (s^)j/ 


(2.8) 


Kq  llgll  2 lIlGgall  2 lllill  2 


Proof.  Given  a,  let  a,  :=  f a.B..  Then,  since 
= li  11 


a,  (x)^  = a a.B.  (x)/(/  x®B.)^^^) 
n 1 1 *'  1 

< £ a^B.  (x)/f  x*^B.  , 

— 11  •'  1 


wc  have 


T t C>  2 2 T 

a Gp3  = ; X dx  £ H a^  = S § • 

Since  Cg  is  symmetric  positive  definite,  this  yields  the  upper  bound  in  (2.7). 


For  the  lower  Ijound,  take  a = (a.),  and  aqain  let  a,  :=  T.  a.B..  By  Lemma  2.8, 

= 1 li  11 


we  have 


(/  x'^B.)-'/^a.  = / x^. 


r.  o(x)a,  (x)dx, 
1 1 1 , P h 

i 


and  lionce 


d k ^ d "4  k. 

(/  x"B.)"^a^  < C(M,Y,k,R)^(t.  - t.)'^t~^  • / x^dx  • / x\  (x)^dx 

X 1 — 1+k  1 iHk  ^ ^ h 

i j 

S+k 

< C(M,Y,k,B)^(t^^j,  - t^)  ^t.^j.  / x*^a^^(x)^dx  . 

*■1 


Thus 
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, . 1 


a.  < C(M,Yik,6)  (t.  - t.)  ^t..^  • / x^B.dx  ’ / x^a  (x)^dx 

X — i+K  1 nk  1 •'  h 

i 


^i+k 

< C(M,y,k,P)^  • k ^ / x®a  (x)^dx  , 

t. 

1 


since  /b,  =(t,  ,-t.)/k.  A summation  on  i 
‘ ' i+k  ’ 


gxves 


a"^a  £ C(M,Y,k,3)^  / x®a^(x)^dx  = C (M , y , k , 8)  ^a'^G^a  , 
which  gives  the  lower  bound  in  (2.8). 

Now,  Deinko  [8]  proved  the  following  result.  If  A is  a band  matrix  satisfying 


-* 

» < 


ii  llyll  q 5.  q .1  llyll  q some  1 £ q £ «>  and  K > 0,  tlien  tlie  entries  of 

A ^ = *“ij*  decay  exponentially  away  from  the  diagonal,  in  the  sense  that  there  exist 
constants  C > 0 and  if  (0,1)  deiiending  only  on  K,  K and  the  bandwidth  of  A 
such  that  l“jjl  £CX^^  Tliis  result  combined  with  (2.8)  allows  us  to  conclude  the 

existence  of  constants  C > 0,  X £ (0,1)  , both  dej'omiing  only  on  M,y,k,f?,  sucli  t)iat 

(2.9) 


(G*^  , . I < 
e xj ' - 


We  can  use  (2.9)  to  estimate  |]d  ^Gg^D||  . Indeed, 

|(d'^g“^D)_|  - I (/  x^B.)■^^^(G■^  _(/  x®B 


B 'ij' 


Li. 


< cx 


l'■■^l(/  x^iV/M  x®B,)-^/2 


— j-fk  jHk  j i+k  l^k  1 


< C • M 


g+k  nk 


Thus 


I|<d"^g'’d)  I £ y xl""3l  + )' 

j=l  ® ^ j=l  j=i£l 


C(M,y,k)'’^^(j/i)’''^/^x'^'^l) 


1/2  -1 

< C’M^'  ((1  - X)  + C(M,Y,k,P)) 


£ 


-8- 


and  hence 


i j=l  0 J 

where  C depends  only  on  M,Y,k,6.  With  (2.7),  this  implies  the  second  conclusion 
of  Theorem  2.1. 

For  the  bound,  we  write  the  matrix  problem  (2.6)  slightly  differently.  L 

A :=  diag((/  B^)  ^ ),  E :=  diag ( ( J Then  we  have 

||Pgu||  = J < 1 aj  / B.  = llAgll  I , 

and 

t 1/2 


|e  ^u||  ^ k [[ull  ^ 


Now,  G^a  = u is  equivalent  to  Aa  = (AG^^E)E  ^u,  so 


(2.10) 


But  Ag”^E  = Ad"^g”^d"^E  - (Ad"^g"^DA~^) (Ad"^E) , so 

Up  p 


-1 


(2.11) 


ll%^:||  ^ < II  (da'^  Gg^DA'^ll^  ||ad"^e||  2 • 


We  can  estimate  each  of  the  latter  two  quantities.  Indeed, 


AD 


^E  = diag((/  x^^B^ 


)^/^  , 


(2.12)  IIad^'^eII  < max  (/  B ) (M,  Y ,)c , B)  t"^.  (/  B.)"^t®_(/  B.) 

2 — , 1 i+)c  ' 1 it)c  ‘ X 


-1  -B 


,-1^6 


;l/2 


l^i^n 

j<  C(M,Y,k,B) 

Now  we  prove  a final  lemma. 


-1 


I,emma  2.7.  Ix!t  A be  an  n x n matrix  wliose  inverse  A ^ = (a  .)  satisfies 


l“ijl  ^ where  c > 0,  0 < X < 1.  Let  D = diaq(d.)  satisfy  d.  > 0 for 


1 < i < n and 


Then  B :=  D 


Proof.  Take 


For  each 


Hence 


whicli  complete 


if  1 < 


d./d.  < 
J 1 - 


K(j/i)  if  1 < 


satisfies  ||  B ||  ^ £ C (c , A , K , m)  . 


X - W<3  have 


= I dT^d.a 

j=l  ^ 3 

Bgll  2 = j (j 

1=1  3=1 

n n 

< I I }'  d.^d^la.  ,|)  • 

“i=l  3=1  ^ 3'  13' 


fixed  i,  we  have 


=1  3 = ] 


< K^c  J + K^i 

j=l 


2 -1  7 

< K c(]  - A)  ^ K c 


C'(c,X,K,m)  . 


IIbxII  I ±C'  (c,X,K,m)  )’  I |u,  . |x^ 

i j ^1  1 

n n 

= C (c , X , K,m)  5^  X , J 
j=l  ^ i=l 

£ C (c,X,K,m)  • c • 2(1  - 
= C(c,X  ,K,m)  l|x||  2 < 

the  jiroof . 


To  apply  tlio  Lemma,  we  note  that 


which  completes  the  proof  of  Theorem  2.1. 


..J'- ... 


S3.  A Ritz-Galerkin  Method 


Consider  the  model  problem,  where  a > l , 


(3.1)  -u"(x)  --u’(x)  = f(x)  on  (0,1)  , 

u(l)  =0,  u'  (0)  =0  . 

Such  an  equation  arises,  with  a = n - 1 , when  one  changes  variables  in  a rota- 
tionally  symmetric  Poisson's  equation  in  b".  Let  us  write  (3.1)  in  the  form 

(3.2)  -(x*^u')'  = x**f  on  (0,1)  , 

u(l)  = 0,  u'  (0)  - 0 . 


1 “ n ? 5 

Let  0 < a < b.  Define  H (a,b)  :=  {v  e L_(a,b)  : / x (v  + (v‘)  )dx  < “), 
~ a 2 ‘ 

°1  1 ^ 

H^(a,b)  :=  H^(a,b)  Cl  {v  « C(a,b)  : v(b)  = O).  We  will  )iave  need  of  the  following 
Poincare-typo  inequality. 

Lemma  3.1.  Let  v f H^(a,b).  Then  , < 2(1)  - a)  ||x"''^v’|]  . 

Proof.  We  have,  integrating  by  parts. 


b b 

t 01,  .2,  2 t ,u^l  a^l 

/ X v(x)  dx  “ - ~ ^ )v(x)v'  (x)dx 


l-cTTT  II  — llL^(a,b)  ll<^“"  -"")’'-"^^-llL,(a,b) 

lair  II>‘“^'''IIl  ,a,b)  - ^“^^»^-“IIl  (a,b)  H — ’ll 

2 oo  2 


,,a+l  II  II  II  a/2  , n 

-a  )b  ||x  v||  ||x 


An  elementary  comjiutation  sliows  (L*^*^^  - a*^^^)b  °'  < (a  + 1 ) (b  - a)  if  a ' 0, 
while  if  a = 0 the  upper  bound  is  trivial.  Thus  the  lemma  holds.  A more  general 


version  of  this  result  may  bo  found  in  (2) . 


For  u,v  f 11^(1),  define  (u,v)  :=  J x'*u'v'dx,  ||u||  :=  (0,11)^^^.  Consider 

Q E _ K h 


the  following  problem. 


Problem  (P)  . Find  u f II  (I)  such  that  for  all  v f H (I),  (u,v)  = / x f(x)v(x)dx. 


-12- 


rttS' 


Let  us  assume  that  x f e L^CI).  It  is  easy  to  show  that  a unique  solution  u 

to  Problem  (P)  exists  and  that  ||uH  ^ 2 l|x°''^^f|| 

E 

We  now  define  the  finite-dimensional  problem.  Let  0 = x^  < x^  < •••  < x^^  = 1 
be  any  partition  of  1,  and  lot  = s[*  be  defined  as  in  §2,  with  1 < v < k - 1, 
Let  :=  {i^  e : i^(l)  = O),  and  consider  the  following  ijjroblem. 

Problem  {P  ) . Find  u e such  that  for  all  v e S^,  (u,  ,v,  )„  = / x°f  (x)  (x)dx . 

n_  h h h h E ^ h 

It  is  again  easy  to  show  a unique  solution  u^  exists  and  is  characterized  by 


(u  - u,  ,v  ) = / x'*  (u  - u,  )'v,'dx  = 0 

h h E •'  h li 


for  all  * S , equivalently 

<3-'’)  \U  - u II  = in^j^  ||u  - V II  , 

v._e  S 
n 

where  u is  the  solution  to  Problem  (P)  . Since  l|v||  = ||) 


jx°'^^v'||  , (3.4)  indicates 

^2 


the  natural  norms  in  tliis  problem  are  the  weighted  norms.  We  will  pursue  this  a 

little  further  and  show  how  Nitsche's  trick  carries  tlirough.  We  need  two  lemmas. 


Y 

Lemma  3.2.  Assume  the  partition  (x^)  satisfies  ^ M(i/i)  for  1 < i j • 

there  exists  a C = C(M,y,a,k)  such  that  for  all  z f il^(J), 

a 

(3.!^)  ^inf,  ||x"/^z  - z^)  II  < Ch  ||x“/\-  II  . 


Proof . Given  z,  we  construct  in  two  stages.  For  the  first  stage  define  a step 

function  t(x)  by  C(x)  = z(x^)  for  x^_j  ^ x^,  1 f i ^ N,  and  CIO)  = C(x^). 

This  makes  sense  because  z f implies  z is  continuous  on  (0,1). 

Now,  for  1 ^ N we  have 

||x°'^^(z  - C)  II  f / x“(z(x)  - z(x  ))^dx  . 

m 

The  integrand  z(x)  - z (x  ) is  a function  in  H^(I  ) which  vanishes  at  x . By 

m am  m 

Lemma  3.1,  the  integral  is  bounded  by  4 (x  - x ,)^  / x”z'(x)^dx.  Summing  from  m = 1 

m n'-l  j 

to  N establishes  the  inequality 

(3.6)  ||x“''^(z  - C)  II  , 1 2h  ||x“^^z‘ II  . 

2 2 


Li 
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The  second  stage  of  the  construction  involves  finding  a smoother  approximation 


z.  to  the  step  function  C.  Let  z.  :=  E.  , C(t.)B..  Then  for  any  x f I, 
h n 1=1  1 1 


C(x)  - z (x)  = [ (c(x)  - C(t.))B.(x)  . 

i=l  ^ ^ 


Suppose  X f . Then,  since  B^  = 0 on  for  i > k,  we  have 


C(x)  - = I (C(x)  - C(t^))B^(x) 


1=1 


= (z(x  ) - z(x  ) )B.  (x) 

i=l 


= 0. 

Now  suppose  m > 1 and  x e I = (x  ,,x  ) = (t,  , , ,,  , ,<t  ,,  , , ).  Then 

m m-1  m (m-l)(k-v)^k  m(k-u)  + l+v 

B.  (x)  ^0  iff  i,  :=  (m  - 1)  (k  - v)  ^ 1 < i < m(k  - v)  + v =:  i,  • Hence 
1 L — — R 


R 

C(x)  - z (x)  = y (C(x)  - C(t.))B.(x)  = y (z(x  ) - C(t,))B,(x) 
h . . 11  . mil 

1=1  


I, 


For  each  i in  the  range  of  summation  wc  have  C(t.)  = z(x.)  for  some  j, 

X ^ ^ 

in 

j ^max  (1  ,in  4 1 - k)  =:  j^.  Hence  “ C(t^)  = / z'(s)ds  (the  integral  exists  because 


X.  > X,  > 0) , so  wc  get 
1-1 


X . 
1 


lz(x  ) - C(t.)|  < (/  s“z*  (s)^ds)^/^  • (/  s~“ 

m 1 — 


ds) 


1/2 


X . 
1 


Hence  for  x r I , 
m 


x°  (C  (x)  - z (x)  ) ^ ^ x°  1 (z  (x  ) - C (t , ) ) ^B.  (x) 


R m 


^ x°  (/  s°z'  (s)  ds)  ■ (/  s °ds)B^(x) 


1=1,  X. 

L 1 


< (|s°^^z'||  ^ , (x/x,  )”(x  -X.  ) ^ B.  (x) 

Ljix  ,x_^)  m 3^  m 3^  J.  1 

m L 

1 L (X,  ,x 

2 

■^m 


"-’II  L,(x,  ,x  • 

23m 
m 
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Integrating  over  I yields 
n 


2 in  2 j m 

m 


Summing  over  m gives 

II  Q/2  \ II  ^ ^ II  I II  2 

II’'  (C  - II  1 - ^ ll’'  ^ II  T ' 

which,  conOainod  with  (3.6),  finishes  the  proof. 

-1 

Lemma  3.3.  Let  f f L^U),  and  define  g via  g(x)  = x / f(t)dt.  Then  g e L^d) 

llgll  L 1 2 II  f II  . 

S ^2 

Proof.  See  [9,  p.  532)  or  (11,  p.  240). 

We  now  show  how  Nitsche's  trick  works  in  our  situation.  Let  u be  the  solution 

of  Problem  (P)  , u.  the  solution  of  Problem  (P._ ) , and  assume  x./x.  < M(i/i)^  for 
h h 3 1“ 

1 £ i i 3 . Then  we  have 

Lemma  3.4.  (u  - u ) ||  f.  Ch  ||  u - u ||  , where  C = C{M,y,a,k). 

n n E 

Proof.  Let  e(x)  :=  u(x)  - u.  (x)  . Let  z(x)  solve  -(x'*?.')'  = x'^e  on  I, 

h 

z(l)  = z'(0)  = 0.  Then 


(z,v)  = / x“o(x)v(x)dx 

^ 0 


for  all  V f For  x c I we  have  the  formula 

-a  a 

(3.8)  z(x)  = j t / a e(s)dsdt  . 


Differentiating  (3.8)  twice  yields 


“/2  „ -1  f a/2  , _ a/2  , , 

X z'(x)  = ax  J (s/x)  s e(s)ds  - x e(x) 
0 


and  hence,  by  Lemma  (3.3),  x'*'^^z"  e Ljd)  and 


|x"/V||  l(2a  t 1)  ||x“/^e|| 

2 2 


Now,  by  (3.3)  and  (3.7)  wc  liave  for  all  z,  c S , 


(3.10) 


I “/2  I,  2 

|x  ell 


(z.e). 


= (z  - z^.e)^ 


1 llell^llz  - z^ll^ 


loll  J|x“/^z-  -z')||^ 


Note  that  the  space  (S  ) ' = (v'  : v e S } is  exactly  equal  to  the  space  S ' 

h h k-l,v 

Furthermore,  it  is  easy  to  sec  from  (3.8)  that  z'  e H^(I).  Thus  we  may  apply  Lemma  3. 

h 

to  find  a function  C.  * S (and  hence  a z.  e S with  z/  = ) such  that 

h K-l,v-l  h h h 


(3.11) 


|x“'^^(z'  - Cj^)  II  £ C(M,Y,a,)<)h||  x“'^^z"||  L^ 


The  relations  (3.9),  (3.10)  and  (3.11)  yield  the  result. 


The  argument  cont.iined  in  I.emma  3.4  obviously  gives  an  error  bound  of  the  form 


|u  - u II  £ Ch  ||x°^^f  II  ; thus  we  have  estimates  for  1|  x 

n b L- 

lx“/2 


■ (u*  - u* ) II  and 
h 


‘ (u  - I • It  is  more  interostiJig  to  consider  the  quantity  ||u  - u jj  for 

n L h L 

2 P 

p = 2 and  p = <“,  as  an  a priori  bound  on  those  would  ensure  the  error  could  not  be 


l>adly  beliavcd  at  x = 0,  at  least  for  reasonably  smooth  u.  Our  first  step  is  to 


bound  u'  - u' . 

h 


Y 

Theorem  3.5.  Assume  the  partition  (x^)  satisfies  xVx^  £ M(j/i)  for  1 ^ i j < 

M > 0,  Y ^ 1 . Then  there  exists  a C = C(M,Y,Q,k)  > 0 such  that 


(3.12) 


1“'  ■ “A"  I - ^ Ih'  - vjl 

h 


and 

(3.13) 


|u-  - ujl  < CmJ/2  inf  ||u-  - v'll  ^ 

00  = ®h  « 


where  M is  the  global  mesh  ratio  defined  after  (2.1). 

= 1 

Proof.  Ix!t  P : L„(I)  -►  S.  . . via  ( x**  (w  - P w)v.  dx  = 0 for  all  v.  f S, 

a 2 k-l,v-l  ^ ah  h 


k-l,v-l 


*h 

Since  (S  ) ' 


S , (3.3)  implies  u'  = P (u').  Thus 

K“J  > V“  1 n Q 
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Theorem  2.1  gives 


for  all  V.  e S . Hence  for  p = 2,“>, 
n 


|u'  - ujl  L ill^  - 11“’  - v,;l 


v^e  S 
n 


- ^ ll^all  p’  11'^’  ■ ''hll  I 

P 

n 


1 1 


< CM^  ^ inf  II  u‘  - vjl  , 

= oh  P 

v^f  S ^ 

n 


and  the  proof  is  complete. 


Corollary  3.6.  (a)  If  u e H'^(I)  where  i £ j 5.  » then 


tu-  -u'll  ICh^-l  ||u<^>||, 

^ ^2  ’'2 


(b)  If  u £ W^'“{I)  :=  {v  : v,v' , . . . e L (I) } where  1 1 j 1 k,  then 


|u'  - u:  II  , < CMy^h^'^  ||u'^’ 


h"  L 


Proof.  The  Corollary  is  an  immediate  consequence  of  (3.12),  (3.13)  and  well-)£nown 


approximation  properties  of  the  spaces  S (sec,  e.g.,  |16)). 

K f V 


We  now  use  an  argument  like  that  of^Lemma  3.4  to  give  an  a priori  bound  on 
|u  - u ■' 


Theorem  3.7.  Under  the  hypotheses  of  Theorem  3.5,  there  exists  a constant 
C = C(M,y,a,k)  > 0 such  that 


(3.14) 


|u  - Uj^ll  < Ch  I|u-  - u-ll 


Proof . Let  e(x)  :=  u(x)  - u^^(x),  and  let  z(x)  solve  -(x  z')'  = e on  (0,1), 


z(l)  = 0.  The  boundary  condition  to  impose  at  0 is  lim  x**z' (x)  = 0.  Then  there  is 

x-K)+ 


a unique  solution  z,  given  by 


I t 

z(x)  = j t '*/  c (s)ds 


(3.15) 


Differentiating  twice  and  using  Lemma  3.3,  we  find  x z « ^2^^^ 
ll>‘°z"ll  T 1 (2a  + 1)  ||e||  . Now, 

h ^2 


/ (x°*z')'e  = j x°z'e'  = / x°(z'  - 


for  all  Zj^  E S , using  integration  by  parts  and  (3.3).  By  Lemina  3.2,  there  is  a 

z f S*’  with  ||x”(z'  - z* ) II  £ Ch  ||x”z"||  , C = C(M,y,n,k).  Thus 

h h 

l|e||  ^ 1 ||x“(z-  - z')  II  ^ lCh||x“z"||  lle'lL 

2 2 2 2 2 

< Ch(2a  + 1)  ||e|l  lie'll  . 

^2  ^2 

Corollary  3.8.  If  u f H^(l)  where  1 £ j ^ then  ||u  “ “u  II  , £ 11“^^'  II  r ' 

^2  ^2 

so  we  have  an  optimal  convergence  rate  in  I.^- 

The  Corollary  is  an  immediate  consequence  of  (3.14)  and  Lemma  3.6. 

We  now  turn  to  the  more  delicate  matter  of  obtaining  an  estimate  for  ||u  - u || 


To  begin,  let  z c (0,1)  and  let  G(x)  - G (x)  be  the  Green's  function  for  (3.1), 

z 


(3.16)  G(x)  = 


1-a  , . 

z -l,0<x<z 


X -l,z<x<l 


(u  1)  , G(x)  = 


log  1/z,  0 < X < z 


log  1/x,  z ^ X ^ 1 


(a  = 1) 


Then  for  all  G,_  c S , 
h 


!(z)  = / x”G'(x)e'(x)dx 
0 


= / x”(G'(x)  - G'  (x)  )e'  (x)dx  , 

0 ^ 


(3.17) 


|e(z)|  < ||e'l|  inf  / x“|g'(x)  - G'(x)|dx  . 

■»  „ °h  0 
G,^€S 
h 


We  have  Theorem  3.5  to  bound  the  first  factor  here,  so  our  goal  is  to  bound  the 


second  factor.  Our  method  will  bo  to  construct  a suitable  "interixjlant"  of  G.  The 


ft 


next  theorem  gives  the  result.  First  let  us  define,  for  a given  partition 
0 « Xq  < < ...  1 k + 1 ^ r ^ N, 


N j ^ 

^ ==  .1^  ’‘jlk  - ’'j-k'  ^j+k-i  ■ 


where  x.  :=  x = 1 for  j > N.  Note  that  if  x.  = (j/N)^,  then 
3 N 3 

Y - k < -1  =>  < C(k,Y)h^(r  - 

Y - k = -1  =>  S < C(k,Y)h^|log(N/(r  - k)  ) | 


Y - k > -1  =>  < C(k,Y)h 

Theorem  3.9.  Let  x,  < z < x,  (1  < J < N)  and  let  G(x)  = G (x)  be  the  Green's 
J“i  u “ 2 

function  given  in  (3.16).  Assume  the  mesh  satisfies  £M(j/i)^,  1 1 i £ j.  Then 

there  exists  a constant  C = C(M,k,a,Y)  > 0 such  that 

1 

(3.19)  inf  / x“|G'(x)  - G*  (x)  |dx  < C(M,k,a,Y)  (h  + S_  . ) . 

®h  0 " 


In  particular,  if  the  mesh  is  quasi -uniform  (i.e.,  x . - x , , > nh  for  1 < i < N, 

3 3-1  - - - 

so  Y = 1),  then 


inf  / x“|g’(x)  - G' (x)  |dx  < C(u,k,a) 

0 " 

G^f  r. 

h 


h log (1/x  ) , k = 2 
u 


, k > 3 . 


Corollary  3.10.  If  u ^ W-^'  (I)  where  1 ^ j < k,  then 


^ C (M,k  ,a,  Y)  M (h^  + )i^  • lo  particular,  if  the  mesli  is  quasi- 


uniform then 


|o||  ^ < C(M,k,a)  ||u' 


irllog  h|  , k = 2 


, k > 3 


Corollary  3.10  follows  immediately  from  Corollary  3.6,  (3.17),  and  Tlieorcm  3.9. 

Its  content  is  tliat  we  have  an  optimal  rate  of  convergence  in  except  for  the  case 

k = 2 of  piecewise  linear  trial  functions  S^.  This  is  the  same  result  as  found  in 
(15).  An  example  will  be  presented  later  to  show  the  | log  h|  factor  cannot  be  removed 


if  k = 2.  First  wo  prove  Tlieorcm  3.9. 
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Proof  of  Theorem  3.q.  The  proof  is  somewhat  lengthy,  but  straightforward.  Wo  construct 
a map  T : such  that  Tf(x)  depends  only  on  the  values  of  f in  a neighborhood 

of  X,  i.e.,  T is  local.  Then  we  let  - TG.  Now,  G consists  of  two  smootli 

pieces  joined  together  continuously  at  k.  Hence  the  problem  of  estimating  G'(x) 
can  be  divided  into  tliree  cases,  depending  on  wliether  x is  sufficiently  far  to  the 
left  of  z,  sufficiently  far  to  the  right  of  z,  or  in  the  vicinity  of  z.  The 
estimates  for  these  cases  are  combined  to  yield  the  result. 

n 

To  begin,  let  ^ L'  witli  X.g  = / f.g  dx,  X.B.  - 6,.  and 

1 1 1 ^ t it]  1] 

_ j i 

ll^-ill  I 5.  - t.)  ; such  functions  f.  are  constructed  in  (41  and  are  used 

A K 1+K  X I 

in  Lemma  2 . 5 above . 

n 

p ) 1 

For  f c 1 define  Tf  = ) (A.f)n..  Then  T maijs  I,  to  S . VVe  would  put 

1 . r 1 '1  ' 

r =1 

°h 

G " TG  excoi  t we  require  G f S . Thus,  define  G = TG,  G = G - (X  G)li  . Because 


h h n n ’ 


B.  (1)  - 0 for  1 i < n,  wc  liave  G.  f S . Now, 
X — h 


1)  ||x“(C'  - G-)|l  < ||x“(G-  - 5;_)||  -t  ||x"(G'  - G'l 

l.et  us  note  for  future  reference  tlio  inequality 


/ |b’  (y)|  dx  < 2 , 


which  follows  from  tlie  equality  (see  (4.f.)  of  111) 

k - 1 _ k - 1 

i>k  t.^.  , - t.  i,k-l  t.,,  - t.  , ”itl,k-l 
rtk-1  I Ilk  nl 

exjiressing  the  derivative  of  tlie  i.''*'  B-spltne  of  order  k in  terms  of  P-:-.i>linos  of 
order  k - 1. 

Wc  can  easily  estimate  tlie  second  integral  on  tlic  ricjlit  liand  side  of  (3.21).  For, 

G - G = (X  G)B  and 
h h n n 


t , 
nik 

|X  g|  = / f (X 

n ^ n 


)G(x)dx  < D,  G(t  ) . 
— k n 


But  t = , > 1 - h and  thus 

n N-1  -- 
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G(t  ) < 
n ~ 


- 1,  a ^ 1 


log  (1/t  ) , a = 1 
n 


< C(a)h  . 


Thus  |A  cj  < C(a,k)h,  and  so 
n 

(3.2 J ) ||x“(g;  - G-)  II  = / x“|X  G|  |b'  (x)  |dx  < 2C(n,k)h  . 

h n ' n II  — 

Now  comes  the  task  of  estimating  tlie  other  term  on  the  right  hand  side  of  (3.21). 

For  fixed  x and  f sufficiently  smooth,  wo  liave  f (s)  = P(s)  + R (s),  where 
k-1  ^ 

r(s)  :=  I (s  - x) ^ f ^ (x)/r!.  It  is  clear  from  the  definition  Tf  = E(A.f)B.  and 

r=0  ^ ^ 

the  fact  that  ~ that  T preserves  polynomials  of  order  k,  hence 

(Tf  - f)'(x)  = (TR  - R )'{x)  = (TR  )'(x)  = [ (X.R  )B!(x)  , 

XX  X r 1 X 1 

1 

so 

(G  - G)  ■ (X)  = [ (x.R^)n;  (X) 

II  .1X1 

1 

where  ts  G minus  its  Taylor  polynomial  expansion  at  tlie  point  x (whicli 

exists  for  x ^ z) . 

Suppose  X e I with  m < .1  - k.  (Recall  z c I .)  Then,  witli 

m ~ J 

i :=  (in  - 1)  (k  - u)  + k and  i : = (m  - 1 ) (k  - u)  + 1 , 


(3.24) 


R 

(G.  - G)  • (x)  - y (X  ,R  )b:  (x)  , 

h . 1 X 1 


since  supii  B.  1 $ if  and  only  if  i < i < i„.  For  i < i.  , 

1 m 1,  — - R — R 

X.R  = / f.(s)R  (s)ds  = 0 since  R (s)  = 0 for  x < z and  s < z.  Thus 

1 X ■'  1 X X 

i 

(G,  - G) ' = 0 on  1 for  m < J - k,  which  is  the  result  for  x sufficiently  far 
11  m ~ 

to  the  left  of  z. 

Now  suppose  X f I with  m > J + k.  Again  we  have  (3.24).  Furthermore,  for 
m “ 


I < 1 < i„, 
L - - R' 
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sup  |r^(s)1  = sup  1/  (x  - i)^  ^G*'^Nc)dC/(k  - Dil 

1/  (X  - C)’'”^d£| 


t ,<s<t. , , 
1 Hk 


1 Ilk 


1 -a-k 

^ C(a,k) t.  sup 

t.<s<t.  , s 
1 Ilk 


< C(n,k)t^  ” ^ max(x  - t,  ,t.  , - x 

- X,  w X,  X +k  m-1 

L L K 

. , , l-o-k  , .k 

' ( (a,k)x  , max(x  - x ,>>:,,,■!<  ,)  < 

— m-k  m m-k  m+k-1  m-1 

,,  , 1-a-k 

3J11C'-'  G (5.)  - c,  i,  and  t.  >x  ,,t.  ,<x 

k,n  — m-k  ^1^'**^  ~ m+k-1 

llenr-f  for  i < x < i , 

I,  --  — K 

f.  (s)R^(s)drl 
^ i 

1>  • C(u,k)x^  ^ ^ niux(x  - X ,x  , , - x ,)*^  , 

k m-k  m m-k  m+k-1  m-1 


and  so 


U.2G) 


n - 

J x'|g'  -G'ldx<  y |.\.R  I / x^jn'(x)|dx 
1 " i ' 1 ' 

m lx  m 


, , rt  1-0 -k  ,k 

< ('(o^Klx  X , mn>:(x  - x , ,x  , , “ x .) 

— m jn-k  m m-k  m-J 

1 - k k 

C(M,k,o,>)x  max(x  - X .1  1 ^ -1  ) 

m-k  m m-k  m-»-k-l  m-1 


whc'J^  V.V  t lir  rior.k  hyi  >ot  l.ri'.  i ;•  x /.x  , < M(m/(m  - k)  ) . 

m in-k  - 

Summing  (^.2«’)  ff'r  .t  + k m N cjivrs 


(3.27) 


/ x"|g'  - cjdx  ■-  C(M,k,n,Y)S^j^l^ 

’'j+k-1 


wtiioli  is  tho  rosult:  for  x suf  f i oionf.l  y far  to  the  right  of 

Finjlly  wi.-  must  consider  ( lio  interm'-diate  case  xf  I ,J-k<m<d  + 

m 

C,  (x)  - G(x)  = )'  U.G)B.(x)  - G(z)  I G(z)  - G(x) 
h +■  1 1 


k.  Write 


= y 1,(G  - G(z))B.(x)  I G(z)  - C(x)  , 


and  thus 


(3.28)  / x“|g^  - G’ |dx  < ^ |A^(G  - G(z))  1 / x“  | B'  (x)  | dx  + / x“lG'(x)|dx 


< y Ia.  (G  - G(z)  ) |x“  • 2 + c h , 
~ ' 1 ' m a 


whore  i , i are  defined  as  before. 
L R 

For  i < i < i„, 

L - - R 


|A.  (G  - G(z)) 


i+k 

“ / 

t. 


(s)  (G(s)  - G(z)  )ds 


and  thus 


= / f . (s)  (G(s)  - G(z))ds 


< D (t  - t.)  z '^jt.  ■ 

— k x+k  1 a ' i+k  ' 


< C(k,a)hz  , 


X\ 

2 )'  |A.  (G  - G(z)  ) lx"  < C(k,a)h(x  /z)  . 

' -I  ' m • ' m 


If  z ^ x^  (i.e.,  .1^2),  wc  have 

(x  /z)"  < (x  /x  < C(M,a,Y) 

m — J+k-1  J-1  — 

and  hence  (3.29)  is  bcjunded  by  C (M,  k ,a , Y)  li  ■ This  coupled  with  (3.28)  proves  (3.19) 


in  the  case  z ^ x^^. 

Finally,  sup)X)se  0 < z < x^^.  In  this  case  define 


G(x)  = 


G (x)  , Xj  1 z ^ 1 

G(x^)  , 0 ^ X < x^ 


By  all  tliat  has  gone  before, 


(3.30) 


inf  / x°1g'  - G' |dx  ^ C(M,k,a,Y)  (h  + S ) . 
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Also, 


(3.31) 


1 , 1 _ 

/ X |g'  - G'|dx<  / x”c  X '^dx  < c h 
0 -Q  “ - “ 


Combininq  (3.30)  .ind  (3.31)  proves  (3.19)  in  tlio  case  0 < z < x^ , and  tlius 
Tlieorem  3.9  is  piovcn.  It  may  be  worthwliilc  to  note  tliat  lettinq  z 0+  proves  tliat 
(3.19)  holds  even  for  z - 0. 

Example.  Here  is  an  example  tliat  shows  t)ie  |log  h|  factor  cannot  bo  removed  from 

the  upper  bound  of  Corollary  3.10  if  );  = 2.  Let  = j/N,  0 £ j ^ N , li  = N \ 

2 

k - 2,  u(x)  = 1 - X . Let  u.  :=  u,  (x  . ) , 0 < j < N . Then  u,  (x)  is  given  by  (3.3), 

2 h 3 - - h 

/ x°  (u  - u,)‘v'dx  = 0 for  all  v t , 

0 

and  u “0.  For  v,  we  may  choose  a function  such  that  v,'  =1  on  1 v'  --  U off 
N h h 3 li 

1..  (Tlio  fact  tlierc  is  only  one  boundary  fiondition  to  impose  on  v,  makes  this 
3 h 

possible.)  Doing  so  we  find 


/ x“-J. 


;i:Li 


dx  - j x'^(-2x)d>;  , 


1 . 
3 


and  hence 


u . - u , , _ , , , ,nt  2 , . , , n + 2 

3 3-]  p .3 ^r_.J_3_r. JJ 

h a + 2 .a+1  , , , , u-i  1 ' 

3 - (3  - 1) 


Since  u ^ = 0,  we  see 


, .rH-2  , . ,nt  2 

= u.  t 2 h^  1 -J3_r..i> 

3“1  3 (1+2  .all  , . , . a-t  1 

3 - ( 3 - 1 ) 


, N ,a  + 2 , . , , n + 2 

2«-L1,,2  j;  3_  -(3-1) 


n ^ 2 


.1  ■“*■1  I ■ 1 \ ^ 

3=mll  3 - (3  - 1) 


In  particular. 


14.  Kxtensions  and  Numoiiral  Kpnulti'. 


The  purpose  of  this  section  is  to  extend  the  results  of  the  previous  section  to 
nonlinear  equations.  Our  results  will  not  be  as  complete  or  satisfying  as  tlior.e  of 
52  and  53;  in  particular,  wo  will  assume  1 < u < 2 and  will  not  obtain  L error 
estimates.  I'inally,  we  will  show  some  numeiical  results. 

No  begin  with  a lemma  which  will  be  useful  later. 

h 


IjOmnia  4 ■ 1 . Kui^poso  x^  ^ c^h  a 


0 where  c^  > 0 is  independent  of  h.  Then 


tliere  exists  a constant  C = C(k,c  ) such  that  for  any  a > 1 and  any  at 

1 ~ K I V 


(4.1) 


2 e 2 

/ x°v'  (x)  dx  ^ C (k  ,c  ) h / a (x)  dx  . 

r\  fx 


h -I  ^ 2 ^ ^ 2 

rrocj^.  ror  any  V t y , | x (x)  dx  ^ since  v'  is  a poiynomaal 


of  order  k on  (0,x^l  and  all  noinr.  art-  otju  ival  cut  cii  a f ini  t.o-dimonsion..l  njaco. 


1 1 1 
f Cl  2 ; u 2 / a 2 

/ xs5  - j xs'  ■*  ] X.  ^ 


^ 2 . a f ^ 2 


a r 2 cx  r 4 

IXj  / V i / w' 


^ Cj  mind  ,r^)h  J v" 


0 


As  a first  step  to  the  iiunlinear  ju-oMrm  wc  investigate  the  effect  ol  lower  order 
terms  on  the  linCi.r  jiroblern.  Consider  the  (-quation 


(4.2) 


-(x”u')'  + x'*q(x)u  - x"f  (x)  on  (0,)) 


u(l)  = 0,  u'  (0)  = 0 . 


a n 

The  equation  (1.2)  cones  from  -u"  - — u'  + qit  = f after  multiplication  by  x . 


We 


"h 


seek  an  a)'iiroximation  tl  t s of  u as  t)io  solution  of 

h 


111  f) 

/ x^U'v'dx  + / x"qtl,  V dx  = / x''"fv  dx  for  all  v « s'  , 
„ h li  „ h li  i h h 


(4.3) 


equivalently, 

(4.4) 


1 1 

/ x'^Cu'  - U')v'dx  + / x°*q(u  - 1).  )v  dx  = 0 for  all 
0 ^ 0 


V.  6 S 
n 


Let  us  make  the  following  assumptions. 

(4.5)  There  arc  c > 0,  Q > 0 such  tliat  (c  - 1)1^  ■'  q„  ' <3(x)  for  all 


0 ■ • ‘ ■ "0  0 ' ''0  - 
2 

X r (0,1),  where  X is  the  smallest  ix>sitive  eigenvalue  of  the  problem 


-(x“v')‘  = X^x“v  on  (0,1),  v(l)  = 0 = v*(0). 

(4.6)  There  are  constants  c^^  > 0 , M > 0 , y ^ 1 such  tliat  x^^  > c^li  and 


x^/x^  < M(j/i)^  for  1 ^ i ^ j. 


By  (4.5),  unique  solutions  u to  (4.2)  and  U.  to  (4.3)  exist. 

Oh  1 

Let  u,  f .5  be  given  by  / x°*(u‘  - u')v'dx  - 0 for  all  v,  e S . Section  3 

h ^ h h li 

gives  us  estimates  on  ||u  - U|^||  in  various  norms,  licnce  to  estimate  ||u  - Uj^||  in 

those  norms  it  suffices  to  bound  ilu,  - U,  ||  . 

" h h" 

The  equations  defining  u,  and  U,  may  be  combined  to  yield 

h h 

1 ^ 1 1 

(4.7)  / x°(u'  - U')v'dx  i j x*^q(u,  - IJ  )v  dx  = / x”q(u.  - u)v  dx 

' h h h „ li  h h „ h )i 


°h 

for  all  V « S . 
h 

°1 

We  also  note  tliat  for  v f II  (I), 

a 


(4.8) 


o '2  ^ 

XV  + 


a 2 


>_l 

1 £ / 


a 

X V 

2 

'"o 

a 

X V 

a ’ 

'2 

^ 2 r 

a 2 

X V 

+ 

(£  - 

1 / 

X V 

0 ^ 

a ' 

' 2 

'2 

X V 

(c  - 

1)  / X 

V 

a ‘2 

X V 


where  c :=  min(c^,l),  and  also 


(4.9)  / x°v  " + / x“qv^  CqXq  / x^v"" 

We  now  state  tlie  main  results  for  tlie  estimation  of  u, 
Theorem  4.1.  Under  the  preceding  assumptions,  we  have 

M ti/2 


(4.10) 


■‘“h-  VllL/li“h-\l'E 


h “h- 


- II  x“^^(u  - U|^ 
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(4.11) 


ll“h  ■ "iJI  11“  • “hll  L^  ' 

(4.12)  l|u  - Uj^ll  < C(c^,c^,Xjj,q,M,a,)c,Y)h(  ||u’  - u,;||  ||u  - Uj^||  . 

Thus  if  we  )iavc  optimal  convergence  rates  for  ||u  - II  ^nd 

lk-u.||^^. 

Proof.  Let  E,  :=  u.  - U.  . Usijig  K.  for  in  (4.7)  gives 

n n n h n 


(4.13)  / x“e^^  + / x“qE^  = / x“q(u^  - u)E^  ll|x“^^E^l|  Q l|x“''^(Uj^  - u 


a/2,.. 


a/2  , 


Hence  (4.10)  follows  from  (4.8)  and  (4.9). 

For  (4.11),  we  have  )jy  Lemma  4.1  and  (4.10), 
'2  ! a'  2 II M 2 


C(k-c^)h“  ! <}  x\^  = llEjl  ^ < C(C„,A^,Q)  ||u  - ujl 


For  (4.12)  wo  use  a Hi  tsche-type  argument.  Let  z be  tlie  (unique)  solution  of 

-(x“z')'  + x“qz  = u - U on  (0,1)  with  z(l)  = 0,  lim  x°z(x)  = 0.  Then 
^ x‘0t 

l!=«^z’’ll  1 C(q>a)  ||u  - U 11  . By  (4.4), 

^2  2 


II  T - u ) 

)i  n h 


= / x“(z'  - z*)  (u'  - O')  t / x”q(z  - z ) (u  - U ) 

•'  h h h h 

°t)  li 

for  all  z,  £ S . By  Lemitia  3 . 2 we  may  clioose  z'  £ S.  , (which  is  oquiv.'^.l  ent  to 

choosing  z f S*^)  sucti  that  ||x”(z'  --  z')||  ^ C (M  ,<> , k , y)  h ||x^z"l|  . We  also  have 

h h I,^  L^ 

llx  (z  - z ) II  i ^ ^ l|x^(z'  - z' ) II  where  A = A(a)  is  tlie  smallest  eigenvalue  of 
h L^  h L^ 

2u  , 2a 

-(x  v')'  - Xx  V on  (0,1),  v(l)  ■-  0 = v'(0).  Hence 

' ^h^l  L ~ ^h^^  L ll^"^"il  ll>^'’^-"ll  L 11“  ■ L 

< C(q,H,a,l.,Y)h  ||u  - UjJ|  ^ ( ||u'  - U'J|  + ||u  - U|^||  ^^)  , 


and  thus  for  sufficiently  small  h 
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/ x'^u'v'dx  = J x^*f  (x  ,u  (x)  ) V (x)  dx 


(4.15) 

'I  'h 

for  all  V e H . Similarly  it  can  be  sfiown  there  i a a unique  u.  f S which  satisfies 
a h 


(4.16) 


/ x^u'v'dx  = / x“f  (x,ir  (x)  ) V (x)dx 
' h h ' li  )i 


for  all  e S . Our  goal  is  to  estimate  u - u,  . 

h li 


L 


8f 


To  do  so  wo  coi  sider  an  auxiliary  iroblcm  cli  fined  as  follows.  Let  q(x)  :=  - — (x,u(x)) 

dU 

and  let  wen  solve 
a 


(4.17) 


-(x'*w’)'  + X^'.IW  - x”(qii  + f(x,u(x))),  w(l)  - w'  (0)  = 0 . 


By  the  assumption  on  f we  have  q ^ (c^  - D'q  and  lienee  (4.17)  has  a unique  solu- 

“h 

tion,  which  is  clearly  u.  Now  consider  the  proMem;  find  u.  f S such  that  for 

h 

all  V.  t S , 
h 

(4.18)  / x“u'v'  + / x^qiiv  = / x"|quv.  t f(x,u)v.  1 . 

hn  nn'  h n 

X 

For  simplicity,  let  us  assume  1 ^ a < 2 anc?  u f H (I).  Then  by  Theorem  4.1, 

(4.19)  ||u  - Uj^ll  ^ ^ + h ||u'  - ujl  1 Ch*'  ||u*''*  H , 

where  C is  independent  of  h.  Notice  that  this  implies 
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(4.20) 


I.  - =J|  ^ ; ch'-  llu"-'  II  ^ 


Let  US  finally  assume  f(x,t)  is  twice  differentiable  with  resi'ect  to  t,  and 
(3^f/3t^j  ^ K.  Tliis  plus  the  previous  assumptions  yield  the  main  result. 


Tlieorem  4.2, 

Let  1 

llu  - ujl 

+ h llu- 

2 

Proof.  For 

all 

n 

t a- 

(4.21) 

/ X u, 

‘^(,11  ] 1 where  C is  independent  of  li. 

°h  . , . . r , 


/ x°Uj^v^  - / x“f(x,u^)Vj^  = / x“[qu  - qu^  + f(x,u)  - f(x,u^)]V| 


h h 


h"  h 


= /x'*(u  - u ) 1-  ~ (x,u)  + — (x,u  -t  0(u 
' h 3u  3u  n 


We  also  have,  for  all  v.  t s 

h 


2 

= - / x*^(u  - u (x,u  + 0(u  - u))v 

‘ h 2 h 1: 

3u 

< K / x“(u  - u )^|v  I . 


(4.22)  / x^u'v'  - / x°f(x,u 

* n n 


= / X^u/v'  - / x“u'v‘  - / x“f(x,u 
‘ h )i  h li  •'  li 


)V, 


/ x“f(x,u  )v 
‘ li  h 


= / - u-)v-  - / X 


"(U. 


(x,u  + 0 (u  - u ) ) 
» n n n 


Putting  - u,  in  (4.22)  yioidy 

n n n 


r a-  - fa-- 

; X u*  (u^  "•  u.  ) ' - / X f(x,u.)(u.  - u.  ) 

h h )»  h h n 


/x"(u;-u.)^- 

/ ■ “l/  fu  ^’''“h  " ®<"h  ■ “h^’ 

> 

/ X (U-  - u^)  + 

(Cq  - 1)^0  ^ x“(Gj^  - u^)2 

> 

c / x“(G;  - u;)2 

« 

where  e :=  min(c  ,1), 

a5T 

in  (4.8).  Thus, 

by  (4.21) 

- u)))v^ 
n 
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and  thus 


C / X (U-  - U')  < / X U-<U^  - Uj^)'  - / X f(x.Gj^)  - u^) 


< K / x“(u  - - uj 


< K Hu  - Uj^ll  J - U|^)ll  ||x“^^(Uj^  - u,,)  II 

- "C  11“  ■ “iJI  L„  II ■ “h>  II  ll>'“^^<“h  • “h’  II  Lj  ' 


||x“/^G-  - u^)  II  < KX'V^  ||u  - GJI  ^ ||x“/2(u  - G^)  II  ^ 

2 "2 

and  also 

ll’‘“^^<“h  ■ “h>  II  L - 11“  ■ “hll  L ll’'“^^<“  - \>  II  L • 

2 “2 

Rccallinq  Leimna  4.1,  we  obtain  with  the  help  of  (4.19)  and  (4.20), 

ll“h  - “hll  ^ ll“h  - “);ll  i C(k,c^).r“/^  ||x“/^G^  - u,^)  II  t ||x“/^G;  - u^)  II  ^^) 

< Ch""^^  llu  - u 11  ^ llx“^^(u  - Uj^)  II 

“ "2 

< Ch-“/V'^-M|u<’^>||J  ■ 


llu  - ujl  + h llu*  - u*  II  < llu  - ujl  t h ||u*  - u;||  t lluj^  - ujl  + ||u*  - uJ 

< Ch’^||u‘’'>|l.  ^ Ch2’^-'-“/^||u‘’‘>|lf 

^2  ^2 

< Ch’^||u‘">||^^  (1  + h'^-^-“/2||Jk)||  ^ 

^2  ^2 

and  the  proof  is  complete,  since  k ^ 2 and  1 5.  “ 1 2. 

Finally,  here  are  some  numerical  results.  The  Ritz-Galerlii  n procedure  was 


programmed  and  tested  on  the  three  nonlinear  problems  which  follow.  The  mesh  was  ta)^en 
to  be  uniform  with  mesh  width  h,  and  was  ta)cen  to  be  cubic  splines 


(k  = 4,  V = 3).  de  Boor's  package  for  calculating  with  B-splines  (3)  was  used  to 
handle  the  spline  manipulations.  The  nonlinear  problem  (4.16)  was  tackled  with  Newton's 


'^■mwpiipipp 


method  and  an  initial  guess  of  0 for  the  solution.  Tlie  integrals  involved  in  the 
linear  problems  to  which  Newton’s  method  leads  were  evaluated  with  a rom)XJSite  Gaussian 
quadrature  rule  using  5 loints  per  interval.  The  iterations  were  judged  to  have 
converged  when  the  residual  was  about  at  the  unit  roundoff  If.'Vel  (whicli  is  roughly  10 
in  double  precision  on  the  Univac  1110).  The  eriors  in  the  h and  1.  norms  were 

2 S' 

then  estimated  by  an  evaluation  at  -10  e<]ually  spaced  points  per  interval.  In  the 

-s 

tables  we  present  the  error  (the  notation  .1714(-!>)  means  .1734  » 10  ),  the  rate, 

which  is  defined  by  log (c (h^)/e (h^) )/log (h^/h^) ) , and  an  estimate  of  the  constant  C 


in  the  inequality 

IIHIp  1 

P 

, Ch^  l|u<^''’| 

1 , defined 

P 

1 by  C _ = 1 

IHIp'’- 

P 

p = 2 and  p = “'. 

For  eacli  cxamiilc  || 

II  was  determined  analytically  and 

II  was  determined 

liy  numerical 

o> 

. integration, 

sj'ecif  ical  ly , 48- 

point  Gaussian 

quadrature. 

Example  1.  a - 1, 

f(.x,u)  ■ 

64  u 
" 40  ■ 

Tlic  solution 

is  u(x) 

= 2 1 og  7 - 

2 log  (8  - >:^) 

Wc  liavo 

II 

11 , T 

||u^"'"||j_  =1.130 

. Eacli  c 

ase  convri  go 

d in  4 iterations. 

'^2 

)i 

Hell  , 

2 

rate 

C 

2 

l|c|| 

<o 

rate 

C 

GU 

1/4 

.1734  (-0) 

.716(-3) 

3388 (-5) 

.70B(-3) 

l/f> 

.374b (-0) 

3.78 

.783(-3) 

7800 (-0) 

3.62 

.895(-3) 

1/0 

.123b(-f.) 

3.00 

.810(-3) 

2805  (-0) 

3 . 50 

.102  (-2) 

1/10 

.5180 (-7) 

3.80 

.830(-3) 

1237  (-0) 

3.07 

.109(-2) 

1/12 

.2530 (-7) 

3.92 

.848(-3) 

0200 (-7) 

3.73 

.llb(-2) 

1/10 

.8174 (-B) 

3.04 

.804 (-3) 

210(.  (-7) 

3.79 

.122(-2) 

1/20 

.3384 (-8) 

3.95 

.874(-3) 

8945(-8) 

3.84 

.127 (-2) 

Example  2.  a = 1,  f(x,u)  - e”. 

This  is  a well-known  exam)ile  (see,  c.g.  |13)).  A solution  (not  unique)  is 
u(x)  = 2 logd  + P)  - 2 logd  •<  Bx^)  , where  B :=  3 - 2>'*2.  The  theoretical  results 
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I 

I 


do  not  include  this  example  ( — has  the  wrong  sign),  but  no  convergence  difficulty 

oU 


was  encountered.  Each  case  required  4 iterations.  We  have  ||u 


.5685, 


h 

IhIL 

^2 

rate 

rate 

C 

to 

1/4 

• 1565(-5) 

.705(-3) 

.3530(-5) 

.128(-2) 

* 

1/6 

.3154 (-6) 

3.95 

• 719{-3) 

.7465 (-6) 

3.83 

.137  (-2) 

1/8 

.1007 (-6) 

3.97 

.726(-3) 

.2403{-6) 

3.94 

.139(-2) 

' 

1/10 

.4150(-7) 

3.97 

.730(-3) 

.9904 (-7) 

3.97 

.140(-2) 

% 

1/12 

.2010(-7) 

3.98 

.733(-3) 

.4791 (-7) 

3.98 

.14K-2) 

f 

f‘ 

1 

Example  3. 

. a 

= 2,  f(x,u)  = 

3/3  + (u  + 

This 

is  a 

modification 

of  an  example  in  113]. 

The  solution  is 

u(x)  = {1 

2 

X 

/3,-^^^  - .'5 

2 

X /2.  Once 

again,  the 

theory  does 

not  cover 

this  exan 

but  no  numerical  difficulty  was  encountered.  Each  case  took  6 iterations.  We  have 

l|u""'|l,  -1. 

2 *» 


ft 

h 

Hell, 

^2 

rate 

S 

Hell 

CD 

rate 

C 

OD 

1/4 

.2107 (-5) 

.636  1-3) 

.52371-5) 

.134  1-2) 

1 

1/6 

.40761-6) 

4.05 

.6231-3) 

.10751-5) 

3.91 

.1391-2) 

1 

f 

1/8 

.12691-6) 

4.06 

.6131-3) 

.34161-6) 

3.99 

.1401-2) 

u 

1/10 

.51531-7) 

4.04 

.6081-3) 

.14001-6) 

4.00 

.1401-2) 

1/12 

.24731-7) 

4.03 

.6051-3) 

.67501-7) 

4.00 

.1401-2) 
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